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We provide a successful fit for proton-neutron scattering below pion production threshold up to LAB energies 
of 350MeV. We use seven high-quality fits based on potentials with different forms as a measure of the system- 
atic uncertainty. We represent the interaction as a sum of delta-shells in configuration space below the 3fm and 
a charge dependent one pion exchange potential above 3fm together with magnetic and vacuum polarization 
effects. Special attention is payed to estimate the errors of the phenomenological interaction. 
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The study of the NN interaction has been a central and re- 
current topic in Nuclear Physics for many years (see e.g. IllQl 
and references therein). The standard approach to constrain 
the interaction is to analyze NN scattering data below pion 
production threshold and to undertake a partial wave analysis 
(PWA), the quality of the fit being given by the ;i^^/dof value. 
Only by the mid 90's was it possible to fit about 4000 selected 
NN scattering data after discarding about further 1000 of 3(T 
inconsistent data with a j^/dof < 1 and incorporating charge 
dependence (CD) for the One Pion Exchange (OPE) potential 
as well as magnetic and vacuum polarization effects yl]. This 
benchmark partial wave analysis (PWA) was carried out us- 
ing an energy dependent potential for the short range part for 
which nuclear structure calculations become hard to formu- 
late. Thus, energy independent high quality potentials were 
subsequently produced with almost identical ;|f^/dof ^ 1 for 
a gradually increasing database iQ-Sl- While any of these po- 
tentials provides individually satisfactory fits to the available 
data it is striking that there are no published error estimates 
of the potential parameters. Moreover it should also be no- 
ticed that the existing high-quality potentials are different in 
their specific form; they range from local to non-local in dif- 
ferent versions of nonlocality. Thus, scattering phase-shifts 
and observable amplitudes are not identical and in fact the 
existing set of high quality potentials as a whole provides a 
distribution of scattering observables accounting for system- 
atic uncertainties in addition to the statistical uncertainties ob- 
tained from the fitted data for each individual potential. Given 
the fact that these interactions are just constrained to the elas- 
tic scattering data (and eventually to the deuteron) which go 
up to the pion production threshold, one is physically prob- 
ing the interaction with a resolution not finer than the shortest 
de Broglie wavelength AX = h / ^Mniu^ ^ 0.5fm. Thus, for 
practical purposes it may be advantageous to consider coarse 
grained interactions. This is actually the physics underlying 
the so-called Viowk approach [8] in which an effective interac- 
tion in a restricted model space is built. By starting from dif- 
ferent high-quality potentials with a common charge depen- 
dent OPE interaction the CM-momenta above A ^ ^JM^niji 
are ehminated by a suitable transformation and a remarkable 



universal interaction is obtained for p < A. Many of the ap- 
plications of such an appealing interaction have recently been 
reviewed |9]. 

On the other hand, when switching from the NN problem 
to the many body nuclear problem the features and the form 
of the interaction are relevant in terms of computational cost 
and feasibility. The lack of knowledge of a precise poten- 
tial form with finer resolution than Ar ^ 0.5fm suggests to 
search for a description of scattering data directly in terms 
of a coarse grained potential sampled at some sensible "thick 
points". Any sampling procedure necessarily redistributes the 
interaction strength and smoothes the potential as compared 
to the zero resolution limit Ar — > implicit in most poten- 
tial approaches and generating the troublesome short distance 
cores. Of course the optimal way to sample the interaction 
is to provide an acceptable J^-fit with the minimal number 
of sampling points; by implementing this minimal sampling 
we just prevent statistical correlations between the strengths 
of the potential at the sampling points. 

In the present work we provide another high-quality poten- 
tial accommodating these desirable features and (unlike the 
previous approaches) we undertake an analysis of the system- 
atic errors. We do this by assuming that the systematic error 
inherent to any specific choice of the potential form corre- 
sponds to individual uncorrelated measurements. Hence we 
may invoke the central limit theorem to undertake the tradi- 
tional statistical treatment to the mean average and the corre- 
sponding standard deviation of Refs. 1 3-7] without any further 
ado. We will use this compilation as our database. 

A convenient representation to sample the the short distance 
component of the NN interaction was already suggested by 
Aviles [llOtl almost 40 years ago in terms of delta-shells which 
for any partial wave (Z, we take as 

V/-^?«-:^L(^«)''5(r-r„) r<r,, (1) 

with jj.„p the reducedpn-mass and to be specified below. 
In the spirit of Refs. 13|-|7[], for r > rc we use the well known 
long-distance tail of the NN potential 

y{r)= VcD (r) + VoPE (r) , r> rc , (2) 
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FIG. 1. np phase shifts in degrees with 7 < 4 as a function of the LAB energy in MeV. The dashed line is calculated wit the fitted potential 
and the points with error bars represent the mean value and standard deviation of the compilation of the PWA 0] and the six high quality 
potentials f^Sf], We fit the energies £lab = 1,5, 10,25,50, 100, 150,200,250,300,350 MeV. 



where Vcd is the charge-dependent (CD), magnetic and vac- 
uum polarization potential of Ref. [5J, and Vqpe is the one- 
pion-exchange potential. 

The solution of the corresponding Schrodinger equation in 
coupled channels is straightforward; for any r„ < r < r„+i 
with < rc we have free solutions and log-derivatives are 
discontinuous at the r„ -points so that one generates an accu- 
mulated S-matrix at any sampling point providing a discrete 



version of Calogero's variable phase equation 1 11]. Although 
this potential is formally local, the fact that we are coarse 
graining the interaction enables to encode efficiently nonlocal- 
ities operating below the finest resolution Ar. Of course, once 
we admit that the interaction below is unknown there is 
no advantage in prolonging the well-known charge-dependent 
OPE tail and other electromagnetic effects for r < r^. The 
low energy expansion of the discrete variable phase equations 
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TABLE I. Fitting delta-shell parameters (A,,)^'^, (in fm ') with their errors for all states in the JS channel and the corresponding ;t;^-value for 
y < 4 in np scattering. We take N = 5 equidistant points with Ar = 0.6fm. — indicates that the corresponding (^n)/'/' = 0. 



Wave 






A3 


X4 


^5 




'So 


2.12(7) 


-0.987(9) 


- 


-0.087(2) 


- 


0.3476 




- 


1.27(4) 


-0.43(1) 


- 


-0.035(2) 


1.1803 


'Pi 


- 


1.23(2) 


- 


0.079(4) 


- 


0.0088 


'Pi 


- 


1.33(2) 


- 


0.053(2) 


- 


0.4371 


'D2 


- 


-0.64(2) 


- 


-0.082(2) 


- 


1.5671 


'D2 


- 


- 


-0.597(8) 


-0.08(1) 


-0.051(4) 


0.6431 


'F, 


- 


2.6(6) 


- 


0.068(7) 


- 


0.1603 


'F, 


- 


- 


- 


0.060(2) 


- 


0.7759 


'G, 


- 


- 


-0.22(2) 


- 


-0.0137(9) 


0.8090 


'G, 


- 


- 


- 


-0.267(3) 


- 


1.6647 


's, 


1.44(3) 


-0.40(1) 


- 


-0.064(3) 


- 




£i 


- 


-1.67(1) 


-0.412(8) 


-0.200(6) 


-0.031(3) 




'D, 


- 




0.52(2) 


- 


0.041(3) 


0.4750 


'P2 


- 


-0.413(6) 


- 


-0.0398(9) 


- 




£2 


— 


0.65(1) 




0.105(2) 


— 




F2 






U. 11(3) 


-0.069(6) 




U. /i/D 






-0.13(4) 




0.010(3) 






£3 






-0.31(3) 


-0.323(7) 






'G, 








0.115(6) 




1.7552 


'F, 




-0.47(6) 




-0.051(1) 






£4 








0.107(3) 














-0.039(6) 




1.4917 



was used in Ref. 1121 to determine threshold parameters in all 
partial waves. The relation to the well-known Nyquist the- 
orem of sampling a signal with a given bandwidth has been 
discussed in Ref. |fl3|. Some of the advantages of using this 
simple potential for Nuclear Structure calculations as well as 
the connection to the Viowk approach have been spelled out 
already Hi I . 



We use the standard LAB energy values usually 
listed in the high-quality potentials, namely E]^ab = 
1,5, 10,25,50, 100, 150,200,250,300,350MeV and fit to the 
mean phase-shift values at those energies with an eiTor equal 
to the standard deviation. In this way we account for the sys- 
tematic eiTors due to the different representations of the poten- 
tials 1 3-7]. We find that they are generally larger than those 
quoted by the original PWA where only statistical uncertain- 
ties where explicitly discussed for a ^jceii potential form I3t]. 
With these data sets and the given energies we undertake a 
phase-shift fit and determine errors using the standard covari- 
ance matrix. 

As expected, we need at most N = 5 sampling points which 
for simplicity are taken to be equidistant with Ar = 0.6fm be- 
tween the origin and = 3fm. This is the minimal number 
which provides an acceptable fit to the data compiled from 
Refs. |l3|-|7|]. Our results for the pn phase-shifts for all par- 



tial waves with total angular momentum 7 < 4 are depicted in 
Fig.[T] The fitting parameters {KiY^i entering the delta-shell 
potentials, Eq. ([Til, are listed in Table |T] with their deduced 
uncertainties. Of course, a definitive assessment on system- 
atic errors would require testing all possible potential forms. 
Thus, the errors will generally be larger than those estimated 
here. We find that correlations among the different (A,„)]f'^/ val- 
ues within a given partial wave channel are unimportant, and 
hence these parameters are essentially independent from each 
other. This is a direct consequence of our strategy to mini- 
mize the number of sampling points. We find that introducing 
more points or equivalently reducing Ar generates unneces- 
sary correlations. Also, lowering the value of r^ below 3fm, 
requires overlapping the short-distance potential, Eq. ([T]), with 
the CD-OPE plus em corrections. 

We determine the deuteron properties by solving the bound 
state problem in the D\ channel using the correspond- 

ing parameters listed in TableU] The predictions are presented 
in table where our quoted eiTors are obtained from propa- 
gating Table U The comparison with experimental values or 
high quality potentials where the binding energies are used as 
an input is satisfactory. This is partly due to the fact that the- 
oretical errors are about 10%. Of course, one may improve 
on this by using the deuteron binding energy as an input as in 
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TABLE II. Deuteron static properties compared with empirical values and high-quality potentials calculations 





Delta shell 


Empirical[I5— 201 


JNijm i |41 


iNijm 11 |4J 


Keidyj |4J 


AV lo pj 


CU-rSonn ^6] 




2.2(2) 


2.224575(9) 


Input 


Input 


Input 


Input 


Input 






fl fl756('5"l 


fl fl?534 


09571 


0.02514 


0.0250 


0.0256 


As(fml/2) 


0.88(3) 


0.8781(44) 


0.8841 


0.8845 


0.8853 


0.8850 


0.8846 


rm{fm) 


1.97(8) 


1.953(3) 


1.9666 


1.9675 


1.9686 


1.967 


1.966 




0.272(9) 


0.2859(3) 


0.2719 


0.2707 


0.2703 


0.270 


0.270 


Pd 


5.7(2) 


5.67(4) 


5.664 


5.635 


5.699 


5.76 


4.85 


(r-l)(fm->) 


0.45(1) 






0.4502 


0.4515 







Refs. M. 

Fitting to phase-shifts to some accuracy needs not provide 
angle dependent scattering ampHtudes to the same accuracy. 
This is often the case when the form of the potential is kept 
fixed, so that the channel by channel fit is usually taken as a 
first step which is afterwards refined by a full fledged anal- 
ysis of differential or scattering observables and polarization 
data [3-7]. To check this issue we proceed as follows. The 
complete on-shell np scattering amplitude contains five inde- 
pendent complex quantities, which we choose for definiteness 
as the Wolfenstein parameters and denote generically as a 10- 
dimensional array (ai, . . . ,010), which could be determined 



directly from experiment as suggested recently 112 111 . We fol- 
low an alternative procedure and construct, out of the high- 
quality analyses, the corresponding mean value of the Wolfen- 
stein parameters,a,(£'LAB, 0), with their corresponding stan- 
dard deviations, Afl, (£'LAB7 0)> for any given LAB energy and 
scattering angle as 



LAB 



0) 



10 

E 

1=1 



Afli(£'| 



LAB 







(3) 



The total 
age over 



X value is obtained as an aver- 
the chosen reference energies /^lab ~ 
1,5, 10,25,50, 100, 150,200,250,300,350MeV and a 
dense 1 "-equidistant sampling of 0-values. The result is 

xVdof=0.78 

which is equivalent to carry a complete J^-fit to the mean av- 
erage scattering amplitude. Again, this is a virtue of the min- 
imal coarse graining of the interactions which adds just the 
needed number of independent parameters. 

To summarize, we have determined a high-quality neutron- 
proton interaction which is based on a few delta-shells for 
the lowest partial waves in addition to charge-dependent elec- 
tromagnetic interactions and one pion exchange and provides 
a good starting point for Nuclear Physics applications. We 
provide error estimates on our fitting parameters account- 
ing both for systematic and statistical uncertainties of present 
day high-quality analyses of neutron-proton scattering data. 
Deuteron properties are compatible with experimental or rec- 
ommended values. Our method allows to coarse grain long- 
range Coulomb interactions in a rather natural way and hence 
to discuss proton-proton scattering data. 
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